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Abstract 

Scenarios with hidden, spontaneously broken, non-abehan gauge groups contain 
a natural dark matter candidate, the hidden vector, whose longevity is due to 
an accidental custodial symmetry in the renormalizable Lagrangian. Nevertheless, 
non-renormalizable dimension six operators break the custodial symmetry and in- 
duce the decay of the dark matter particle at cosmological times. We discuss in 
this paper the cosmic ray signatures of this scenario and we show that the decay 
of hidden vector dark matter particles generically produce an intense gamma ray 
line which could be observed by the Fermi-LAT experiment, if the scale of custo- 
dial symmetry breaking is close to the Grand Unification scale. This gamma line 
proceeds directly from a tree level dark matter 2-body decay in association with a 
Higgs boson. Within this model we also perform a determination of the relic den- 
sity constraints taking into account the dark matter annihilation processes with 
one dark matter particle in the final state. The corresponding direct detection rates 
can be easily of order the current experimental sensitivities. 



1 Introduction 



One of the most striking features of the dark matter (DM) particle is its longevity 
at cosmological scales. This fact could be accommodated ad hoc by imposing a new 
symmetry (discrete or continuous) which prevents the decay of the dark matter particle, 
such as the i?-parity in the Minimal Supersymmetric Standard Model or the Z2 symmetry 
assumed in many phenomenological models. The longevity of the dark matter particle is 
more elegantly explained, however, if it arises as the result of an accidental symmetry of 
the Lagrangian, in complete analogy to the longevity of the proton, which in the Standard 
Model framework is explained by being the proton the lightest particle carrying baryon 
number. A simple implementation of this idea consists in extending the Standard Model 
gauge group with a non-abelian gauge symmetry, under which all the Standard model 
particles are singlets, which is spontaneously broken by the vacuum expectation value of 
a standard model singlet scalar particle. The renormalizable part of the Lagrangian of 
this model presents an accidental custodial symmetry which prevents the decay of the 
hidden vector bosons, thus predicting the existence of a particle with the correct dark 
matter properties [U |2] . 

The simplest example of such class of models introduces in the hidden sector an extra 
S'?7(2)hs gauge group plus a complex scalar doublet of this gauge symmetry, which ac- 
quires a vacuum expectation value. After the SU (2)hs spontaneous symmetry breaking, 
the renormalizable part of the Lagrangian presents a 5*0(3) custodial symmetry which 
makes the three components of the 5'[/(2)hs vector boson degenerate in mass and stable. 
For wide ranges of the parameters of the model, the relic abundance of the vector multi- 
plets can reproduce the observed dark matter abundance. Furthermore, these parameters 
are also consistent with the present constraints from electroweak precision measurements 
and from direct dark matter searches. 

Being the 5*0(3) custodial symmetry an accidental symmetry it is plausibly broken 
explicitly by higher dimensional operators in the Lagrangian. This is again in complete 
analogy with the baryon number violating dimension six operators that necessarily ap- 
pear in the Standard Model Lagrangian, unless the baryon number conservation arises 
as a residual symmetry of an underlying gauge group. Indeed, there are dimension six 
operators which violate the custodial symmetry which can induce the decay of the dark 
matter particle, whereas analogous dimension five operators are absentj^ The scale of 
custodial symmetry breaking has a lower bound stemming from the requirement that 
the dark matter lifetime has to be longer than the age of the Universe, tdm ^ 10^^ s. 

^In this respect, to consider an accidental symmetry is different from justifying the dark matter 
stabihty from a residual discrete subgroup of a Grand Unification gauge group [21 21 |5J IHl [3 [S] , where 
the latter forbids any decay. 
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Furthermore, the dark matter decay produces a flux of stable particles, such as positrons, 
antiprotons and gamma rays. The requirement that the exotic flux of cosmic rays does 
not exceed the measured fluxes translates into a more stringent constraint on the dark 
matter lifetime and in turn on the scale of custodial symmetry breaking. 

On the other hand, a series of experiments measuring high-energy charged cosmic 
rays have recently reported strong indications for the existence of an excess of positrons 
at high energies. Namely, the PAMELA measurements of the positron fraction show an 
energy spectrum which rises steeply at energies 7-100 GeV, possibly extending towards 
higher energies [9], while the secondary positron flux calculated from state-of-the-art 
propagation models [lO], together with the total electron plus positron flux measured 
by the Fermi collaboration [11], predict a positron fraction which decreases monotoni- 
cally with the energy. Furthermore, the electron plus positron flux measured by Fermi is 
harder than expected from conventional diffusive models [12], also suggesting the exis- 
tence of an excess of electrons and positrons at higher energies, with a cut-off at around 
1 TeV as observed by the H.E.S.S. collaboration [131 E]- It is interesting to mention 
that one possible explanation for the electron/positron excesses is precisely the decay of 
dark matter particles in the Milky Way halo with a mass in the TeV range and a life- 
time ~ 10^^ s [151116]. Nevertheless, irrespectively of the origin of the electron/positron 
excesses, these measurements set a constraint on the exotic flux of electrons/positrons. 
Furthermore, models of dark matter decay are severely constrained by the PAMELA 
measurements of the antiproton-to-proton ratio pjT], which does not show a deviation 
from the expectations of conventional production mechanism. Interestingly, for TeV mass 
dark matter particles which decay via dimension six operators, this value for the life- 
time naturally arises if the dimension six operators are suppressed by a mass scale close 
to the scale of grand uniflcation [181 UHl I2Q1 ISH [22], |23] . Therefore, present cosmic ray 
observations provide a way of testing some scenarios of Grand Uniflcation. 

One of the well-known smoking gun signals for dark matter in the sky is the possible 
observation of a sharp gamma-ray line f2M 123 I2S]- The 77 and jZ lines have been 
exhaustively studied in supersymmetric models, for the neutralino dark matter [271 EH 
[291 [30[ [3T| 132] . in the inert doublet model [33] and for Kaluza-Klein dark matter [M] . 
Recently it has been pointed out that a monochromatic gamma line can be also produced 
accompanied to an Higgs boson [3S]. Note that the gamma lines in all these models are 
induced by annihilation of dark matter particles at one loop level (for an exception see 
Ref. [36]). However, in many situations the disentanglement of the gamma lines from the 
background requires either a dark matter mass in the TeV range or astrophysical boost 
factors to make the signal strong enough. On the other hand, intense gamma-ray lines 
can also appear at tree level in the decay of dark matter gravitinos in supersymmetric 
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scenarios where R-parity conservation is not imposed [371 EHl EH HQl HI] . In both cases, 
the observation of a gamma-ray hne would be a very clean indirect detection signature 
for annihilating or decaying dark matter and is a promising signature to search for. 

In this paper we will analyze the cosmic ray signatures stemming from the decay 
of a hidden SU (2) vector boson. A decay, unlike most annihilations processes, does not 
require any "boost factor" in order to lead to cosmic ray rates well above backgrounds. 
Moreover it leads to larger extragalactic fluxes. We will work out the possible decay 
modes in detail and concentrate on gamma-ray line and anti-matter signatures. Most 
interestingly, we find that hidden vector dark matter decay modes with gamma-ray lines 
in the final state are automatically present already at tree level. These decay modes exist 
for every possible dimension six operator leading to the hidden vector dark matter decay 
and are hence a robust prediction of the model. Taking advantage of the fact that the 
dark matter has spin-1, the gamma lines arise directly at tree level from DM decays to 7/1 
or 'jT] where h and 77 are the standard model and hidden sector Higgs boson respectively. 
Furthermore, one of the operators generates a kinetic mixing between hidden sector and 
the hypercharge f/(l)y giving rise to two-body decay modes into charged leptons. We 
will discuss these properties in light of recent and future cosmic-ray observations. 

The paper is organized as follows: in section 2 we will introduce the hidden vector 
dark matter model and discuss the possible decay modes in detail. In section 3 we will 
discuss these decay modes for the different operators separately and show results for 
several benchmark models, giving particular emphasis to gamma-ray line signatures. In 
section 4 we improve the calculation of the relic density taking into account the dark 
matter annihilations with one DM particle in the final state, and compute the direct 
detection cross sections it gives. We draw conclusions in section 5. 

2 Hidden SU{2) model with custodial symmetry 
breaking 

We consider an extension of the Standard Model where the gauge group contains a hidden 
non-abelian group, S'?7(2)hs, with gauge bosons A^. We assume that this symmetry is 
spontaneously broken via the vacuum expectation value of a complex S'[/(2)hs doublet 
scalar field, 0. We further assume that all the Standard Model particles are singlets 
under S'f/(2)Hs, thus the Standard Model only couples to the hidden sector via the 
Higgs portal term being H the Standard Model Higgs doublet (note that the 

kinetic mixing of the 5'[/(2)hs gauge multiplet with the SM gauge bosons is forbidden 
by the non-abelian character of the extra gauge symmetry). Under these assumptions. 
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the renormalizable part of the Lagrangian reads: 

£ = C'^' - ^F'^" ■ F^, + {V^<P)\V^<P) - \^<P^H^H - ^J0t0 _ A^(0t0)2 , (2.1) 

where = d^(j)—i^T-A^, being r", a = 1, 2, 3 the generators of the hidden SU{2) gauge 
group. If yLi^ < 0, the hidden sector scalar field acquires a vacuum expectation value, v^, 
and the S'?7(2)hs symmetry is broken spontaneously, with = (— /i|A + (AA<^ — 
A^/4). In the unitary SU{2)iis gauge the Lagrangian of the theory is: 

C = CsM-\F,,-F'^'^ + ^{g^v^rA,.A^ + ^glA,.A>^i'+^-glv^^ 

+ lidXr -^iv' + v.fH^H - ^(V + v,f -^{^' + v,r , (2.2) 

which gives Ma = g^V(j)/2 and where rj' is the hidden sector Higgs boson. This Lagrangian 
has only 4 independent parameters, which can be taken as g^, v^f,, A^ and A^- 

Once the electroweak sector is broken, the hidden sector 77' mixes with the standard 
model Higgs boson h' through the Higgs portal interaction A^ 

h' = cos (3 h + sin (3 r] , (2 3) 

Tj' = —sin (3 h + cos (3 r] . 

The complete Lagrangian in the h, t] physical state basis can be found in Ref. [1] as a 
function of g^, f^, A<^ and A^, together with the corresponding expression for the mixing 
angle (3. 



The Lagrangian in Eq. (2.2) has a remarkable property: it displays a S0{3) custo- 
dial symmetry in the component space, which prevents any decay to 5*0(3) singlets 
(such as Standard Model particles or r]'). Consequently, if the model is described just 
by the renormalizable Lagrangian, the three A^ components are degenerate in mass and 
are absolutely stable. Nevertheless, since this 5*0(3) global symmetry is accidental, one 
expects in the Lagrangian the existence of non-renormalizable operators suppressed by 
a large scale A which break the custodial symmetry. The complete list of operators with 
dimension smaller or equal than six which lead, after the spontaneous symmetry break- 
ing of 5*f/(2)Hs and SU{2)l x f/(l)y, to the breaking of the 5*0(3) custodial symmetry 
reads: 

(A) ^V^<P^<PV,H^H, (2.4) 

(B) j^V,<P^^H^V,H, (2.5) 

(C) V^<j>^V,<p F^^^ , (2.6) 

(D) (2.7) 
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In turn, the breaking of the custodial symmetry leads to the decay of the dark matter 
hidden gauge bosons. Let us discuss for each case the dominant decay modes: 



Case A. After the spontaneous breaking of the gauge symmetries the non- 
renormalizable part of the Lagrangian has two parts: = + Cj^. The first 

one reads: 

^A? = ^(^4W + 2Vt;, + ^a^(/.'a,/^' + ^a,/^')) +h.c., (2.8) 



which induces the decay of the dark matter particle into scalars, by means of Eqs. (2.3): 
A — )■ rji], A ^ hi] and A — )■ hh. The corresponding decay rates are: 

11 0^ 2 J {Ml - 4M2)3 

^ 3167r256An ^ 'P J 



r(A^M = 3^^(^^in2/3 + 4t;^,cos^/3 

^(--) - ^iiJ.("^^^"----^)^^^^¥^- (-) 

with A(M^, mi, 1112) = M\ + ml + m\ — 2{m\ + m'^)M\ — 2mlml. 

In addition, there is a second term in the non-renormalizable Lagrangian: 

^If = ^ (^4(^V + 2v\ + vl) '-^B,{h'h' + 2vh' + v'))+ h.c. , (2.10) 

which induces decays into a gauge boson and the hidden sector and standard model 
Higgs bosons, A — )■ 777, A — )■ Zi], A ^ 'yh and A — )■ Zh, with rates: 



T{A ^ 777) 



3 16n 64A4 



11 e^glsm^ew{cosPv^v^ + sm^vvlf 

T(A Zri) = — — X 

^ 3 647r 64A4 



{Ml - M2 + Ml)\ ^X{Ma,Mz,M, 



MlMl J Ml 



(2.11) 



The decay rates into 7/i and into Zh are obtained using Eqs. (2.3) and Mr, — > Mh 



Case B. At low energies the non-renormalizable Lagrangian for case B is = — 



where C^i ^^^^ given in Eqs. (2.8) and (2.10) (note the change of sign in 



the operator involving two gauge bosons compared to the case A). Thus, the decay modes 
and rates are identical to the case A. 
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Case C. After the symmetry breaking the Lagrangian for this dimension six operator 
reads as 

= j^2Al{d^B, - d,B,)d^r]'{v' + v<l>) + h-c. . (2.12) 
Remarkably it leads only to the following two body decays 

1 1 QI^I COS^ 3 COS^ 6w / n no a a 1 

^f-^-^"' = 3ISt4i {(Ml-M^-f^^JMUM^)(Ml-Ml 

r(A^Z,) = i _L S>1 '^00' -'f I ^ 
^ " 316ir KA'Ml V > ^' " 

X [ml - APJ - 3MUMI - 2A/^) - MMMl_l3ll'j . (2.13) 

together with the corresponding decay channels into 7/1 and Zh, which can be derived 
using Eqs. (2.3) and substituting — )■ Mh- 



Case D. The non-renormalizable part of the Lagrangian contains two terms: = 
+ jC^. It contains a kinetic mixing term: 

Cl^=^AF'^'^\ (2.14) 

which leads to decays into fermion pairs, into W~^W~ , into Zt] and into Zh. 

For the two-body decay into fermion pairs we obtain (neglecting fermion masses) 

m ^ff) = l^^-^ ((4)' + (4)') Ma , (2.15) 

where we have introduced a color factor which is C = 1 for leptons and C = 3 for 
quarks. The effective vector and axial couplings to neutrinos, charged leptons, up-type 
and down-type quarks are given by 

d-v = d\ = —d-A = d\ = —d^ = ~2jif2 _ ]\^2 ' (2.16) 
1 4 . 2, \ A/| 5 



w /'2 1\ M2 1 

4 = (2.19) 



The other decay widths are given by 



16A4 12 ^'^"^yMwJ \Ml-Ml 
Ml Ml\ ( 4M2 



z 



and 



T{A^ Zf]) 



1 1 C/2]Vf 2 gjj^2 p 



3 2567r 
X TlO + 



cos^ 6 
M2 



Ml 



— sin 9w 
2/ Aj-2 n]\/r2 nn/T'i 



w K^Ma\MI-M\ 



Ml 



A.Ml. M^M2 - 2M1 - 2M| 



Ml 



X 



X 



X 



1 _ (M^ + M,)VM1)(1 - {Mz - M.y/M^ 



{2.21] 



For the corresponding decay channels into Zh, the last equation holds with the 
replacement of the physical h boson from Eqs. (2.3). 



The second term, on the other hand, reads: 

= j^v^d,Alid,B, - d,B,W , (2.22) 

which leads to two body decays involving the hidden sector Higgs r], using Eqs. (2.3), 
and a gauge boson, with rates: 

m^,,) ^ 3^xtcos^^.cos^/5^^^, (2.23) 
T{A^Zr]) = lJ-^sin2^;^cos2/3x (2.24) 



X ((Ml - Ml + Mlf + 2MlMl) ^ ^ ^3^' ''^ . 



The decay into Zh and 7/1 can be obtained with the substitution cos /3 — )■ sin (3. 

If /3 7^ and /3 7^ 7r/2, then we have interference between the decay coming from 
and C^^, which lead to the corrections: 

I ?> vlg sin OwMz sin 2/3 / M| 



5V{A^Z^) = — - — 

^{Ml-Ml^Mi)^M^™. (2.25) 



^v(A^7h\ ^ 3 t^^ ^singH/Mzsin2/3 ^ M| • . \ 

^ ' 3 647rA4 cos^vK \Ml - Ml ) 

x(Ml-M,;.M|)«™I. (2.26) 



3 Cosmic ray signatures of hidden vector dark mat- 
ter decay 

In this section we firstly give a short introduction to the propagation of gamma rays and 
charged cosmic rays through the Galaxy, and secondly discuss the typical cosmic-ray 
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signatures of hidden vector dark matter, including strong gamma-ray lines and possible 
contributions to the anti-matter fluxes. 



3.1 Cosmic ray propagation 

The hidden gauge boson decay produces a high energy flux of stable particles, such 
as gamma rays, electrons, positrons, antiprotons, neutrinos and antideuterons. The 
flux of high-energy cosmic rays depends essentially on the scale of custodial symmetry 
breaking, which is thus constrained by the requirement that the predicted fluxes do 
not exceed the observed fluxes. We will show that a typical signature of the hidden 
vector dark matter model is a prominent gamma-ray line, with values of the custodial 
symmetry breaking scale close to the Grand Uniflcation scale and in reach of the Fermi 
LAT dark matter searches. Furthermore, as mentioned in the introduction, experiments 
measuring the positron fraction and the total electron plus positron flux indicate the 
existence of an additional source of electrons and positrons at high energies, but no 
additional source of antiprotons. We will also explore the possibility that the decay of 
dark matter hidden gauge bosons could be the origin of such excesses. 



The production rate of particle i = e,7,j9 per unit energy and unit volume at a 
position r with respect to the center of the Milky Way is given by 

Q.(E,n-j/^'^, (3.1) 

where dNi/dE is the corresponding energy spectrum of particle i produced in the decay, 
and p{f) is the density proflle of dark matter particles in the Milky Way halo. We will 
adopt in this paper the spherically symmetric NFW density proflle [12] for deflniteness, 

Pir) oc ui ! / / M2 ' (3-2) 

normalized to 0.3 GeV/ cm^ at the position of the sun, r = 8.5 kpc, although our results 
do not depend much on the speciflc form of the halo proflle. 



Gamma-rays, contrary to electrons, positrons and antiprotons, which will be dis- 
cussed below, do not diffuse and carry information about their spatial origin. The gamma- 
ray signal from dark matter decay consists of several components. The most important 
one is related to the prompt radiation {e.g. flnal state radiation) produced in the decay 
of DM particles inside the Milky Way halo. It depends on the halo density proflle, and 
although the halo proflle is expected to be approximately isotropic, the corresponding 
flux at Earth exhibits a dipole-like anisotropy which is due to the off-set between the Sun 
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and the Galactic center and which can be as large as 20-30% for dark matter lifetimes 
of the order of lO^^s In contrast, the extragalactic prompt component of the 7-ray 
signal, which stems from the decay of dark matter particles at cosmological distances, is 
isotropic. At energies around 10 GeV or below, the magnitude of the halo and extragalac- 
tic fluxes are comparable when looking in direction of the anti-galactic center, whereas 
at higher energies around and above 100 GeV the inelastic scattering between 7-rays and 
the intergalactic background light reduces the extragalactic component considerably (see 
Ref. |13] for a discussion). In the plots we assume a 10% energy resolution, and we show 
also the H.E.S.S [131 IE] results for the electron -|- positron (+ gamma) flux in the TeV 
region, which acts like an upper bound on the isotropic flux. 

For details about the adopted calculation of the electron, anti-matter and gamma-ray 
fluxes we refer to Refs. [151 1^ . 

Electrons and positrons from dark matter decay loose their energy mainly via inter- 
action with the Galactic magnetic field and the interstellar radiation field (ISRF). In the 
first case (assuming injection energies of the order of 1 TeV) synchrotron radiation in the 
radio band with frequencies O{0.1 — 100 GHz) is produced and potentially observable 
(see e.g. Refs. [HI HS], |16] ) . In the second case, the inverse Compton scattering (ICS) of 
electrons and positrons on the ISRF (which includes the cosmic microwave background, 
thermal dust radiation and starlight) produces a second component of gamma rays with 
energies between 100 MeV and 1 TeV [ITJ HE] . However, since electron and positron 
fluxes are always relatively weak in the decay channels we consider we will neglect ICS 
radiation throughout this work for simplicity. 

After being produced in the decay of dark matter particles, electrons and positrons 
scatter on irregularities of the Galactic magnetic field before reaching the Earth, which 
results in a wash-out of directional information. Their propagation is commonly described 
by a diffusion model, whose free parameters are tuned to reproduce the observed cosmic- 
ray nuclei fluxes [19]. As propagation parameters we will adopt the ones of the MED 
propagation model defined in [50], which provide the best fit to the Boron-to-Carbon 
(B/C) ratio: 6 = 0.70, Kq = 0.0112 kpc^/Myr and L = 4kpc. Our conclusions, how- 
ever, are rather insensitive to the choice of propagation parameters. The astrophysical 
background in the e^-channel is mainly due to primary electrons, which are presumably 
produced in supernova remnants, and due to secondary positrons, produced in the in- 
teraction of cosmic-rays with the galactic gas. For these background fluxes we adopt the 
"Model 0" presented by the Fermi collaboration in [12], which fits well the low-energy 
data points of the total electron plus positron and the positron fraction. We allow, how- 
ever, for a 10% rescaling of the electron background in order to improve the agreement 
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of the total flux to the data. 

The antiproton propagation in the Galaxy is analogous to the propagation of elec- 
trons and positrons. However, since antiprotons are much heavier than electrons and 
positrons, energy losses are negligible. Furthermore, antiproton propagation is affected 
by convection, which accounts for the drift of antiprotons away from the disk induced by 
the Milky Way's Galactic wind. For predictions of the antiproton flux we show an error 
band, corresponding to the MIN and MAX model of Ref. [50]. In our plots we present 
actually the p/p-raXio, where we adopt the proton and anti-proton backgrounds from 
Ref. [51]. 

For both, electrons/positrons and anti-protons, the fluxes at the top of the atmo- 
sphere can differ considerably from the interstellar fluxes at energies smaller than ~ 10 
GeV, due to solar modulation effects. To take this effect into account, we adopt the 
force field approximation [52] with (pp = 550 MV |53j . 

The main background in the 7-ray channel is the diffuse emission of our Galaxy, which 
is mainly due to interactions of cosmic rays with the galactic gas and the ISRF |54j . 
This component is by far strongest in the galactic disk region, and it turns out that 
exotic fluxes from dark matter decay would dominantly show up at higher latitudes, 
away from the disk. For this reason they could be misidentified as contribution to the 
extragalactic gamma-ray background (although they can be distinguished by their large 
scale anisotropy, see Ref. [13]). In this work we will show predictions for the averaged 
gamma- ray flux in the region < / < 360°, 10° < |&| < 90°, which offers the best strategy 
for searching gamma-ray lines from dark matter decay [55] • 



3.2 Gamma-ray lines 

The existence of two-body decay modes with gamma-ray lines in the final state are a 
generic prediction of hidden vector dark matter models. We will discuss this for each 



possible operator separately, Eqs. (2.4)-(2.7). 



Case A and B. In cases A and B, Eqs. (2.4) and (2.5), the dark matter particle decays 



either into two scalar particles [rji], hrj, hh) or into a gauge boson and a scalar particle 
(777, jh, Zr], Zh). Whether the dark matter particle decays preferentially into two scalar 
particles or into a gauge boson and a scalar particle depends on the details of the model. 
In both cases, the fragmentation and decay of the Higgs boson or the hidden sector rj 
boson could produce a sizable flux of electrons, positrons and antiprotons. Unfortunately, 
the electrons and positrons produced in fragmentations cannot explain the PAMELA and 
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Benchmark 


Ma 


9ct> 




Mr, 


Mh 


sin/3 


1 


300 GeV 


0.55 


1090 GeV 


30 GeV 


150 GeV 


^ 


2 


600 GeV 


0.6 


2000 GeV 


30 GeV 


120 GeV 


^ 


3 


14 TeV 


12 


2333 GeV 


500 GeV 


145 GeV 


^ 


4 


1550 GeV 


2.1 


1457 GeV 


1245 GeV 


153 GeV 


0.25 



Table 1: Benchmark points used for the calculation of cosmic-ray signatures. 



Fermi excesses and moreover the PAMELA measurements on the antiproton-to-proton 
ratio set very stringent constraints on possible new sources of antiprotons. Interestingly, 
even if the scale A is increased in order to avoid an antiproton excess, the generically 
present gamma-ray lines can still be intense enough to be observed in experiments, due 
to the enormous sensitivity of dark matter line searches. 

This is illustrated in Figs. [T] and |2| where we show the predictions for the positron 
fraction, total electron plus positron flux, antiproton-to-proton fraction and gamma- 
ray flux for two generic scenarios, namely the benchmark points 1 and 2 defined in 
Tab. [Tj These choices of parameters can successfully reproduce the observed dark matter 
abundance and are consistent with all present laboratory constraints. We also show in the 
plots for the positron fraction the results from PAMELA [9J, HEAT [56j, CAPRICE [SZ] 
and AMS-01 |SB]; for the total electron plus positron flux, the results from Fermi |llj . 
H.E.S.S. [ini E], BETS, PPB-BETS [59], ATIC [60], HEAT, CAPRICE and AMS- 
01; for the antiproton-to-proton ratio, the results from PAMELA ^7], BESS95 [61], 
BESS95/97 [62], CAPRICE94 [63j, CAPRICE98 [61] and IMAX [65] and for the gamma- 
ray flux, the preliminary data from the Fermi-LAT in the region between 10° and 90°, 
as well as the extraction of the extragalactic flux from these data [66] . In the gamma-ray 
plot, we also show the H.E.S.S. results for the electron + positron (+gamma) flux at 
high energies, which yields also an upper bound on the overall isotropic gamma-ray flux. 

The branching ratios for these two benchmark points are listed in Tab. |2} Benchmark 
point 1 is characterized by large branching ratios into gauge boson and Higgs, being the 
decay into a monoenergetic gamma line the dominant channel. On the other hand, since 
kinematically allowed, benchmark point 2 is characterized by a large branching ratio into 
two Higgs bosons. It is interesting that, even though the decay mode into monoenergetic 
gamma rays is subdominant in this benchmark point, the gamma-ray line still is a very 
prominent feature in the gamma-ray energy spectrum. 

We estimate that, in the limit ^ f , /3 — )■ 0, the decay rate into 'jh is given by: 

12 



Benchmark 


rjrj hrj 


hh 


IV 


Zrj 




Zh 


1 


- 0.09 




0.04 


0.02 


0.65 


0.20 


2 


- 0.04 


0.62 


0.002 


0.003 


0.15 


0.18 


3 


- 0.04 


0.80 


3 X 10-6 


0.002 


0.0003 


0.16 



Table 2: Branching Ratios for Case A. 



The Fermi-LAT observations of the region |6| > 10° plus a 20° x 20° square around 
the Galactic center constrain the dark matter lifetime to be longer than a few times 
10^^ s at energies below 200 GeV [67] , which is taken into account in the results shown 
for benchmark point 1 (see Fig. [T|, where the line is around 110 GeV. Thus present 
experiments can probe values of the scale of custodial symmetry breaking close to the 
Grand Unification scale. In case of benchmark point 2 the line occurs at an energy scale 
above the ones probed by Fermi, so that smaller lifetime are allowed experimentally. We 
show results for a lifetime 1.1 x 10^'' s, where the contributions to the diffuse gamma-rays 
around 10 GeV and the anti-proton fluxes can be sizeable. The gamma line in this case 
is huge and should be seen by any experiment sensitive to these energies. 



Case C. This operator, see Eq. (2.6), predicts decays only into a gauge boson and a 
scalar particle, either h or rj. A large decay branching ratio into monoenergetic gamma- 
rays is predicted unavoidably, as illustrated in Tab. |3] for the four different benchmark 
scenarios. In the limit <^ M^, the decay rate into 717 is given by 

r,^ ^ ^ 2.7 X 10-s (,,^V . (3,4) 



^4 X 1015 GeV J V Ma J 
which can make the gamma-ray line observable at the Fermi-LAT for values of the scale of 
custodial symmetry close to the Grand Unification Scale, especially for large dark matter 
masses. The cosmic ray signatures of benchmark point 1 for case C are very similar to 
case A, cf. Fig. [!} On the other hand, we show in Fig. [3] the predictions for benchmark 
point 3 with a very large dark matter mass of 14 TeV, which predicts a strong line at 
very high energies and only small contributions to positrons and anti-protons. 

One feature of the model that is in principle present for each operator, and which 
we want to illustrate for case C, is the general existence of two independent gamma-ray 
lines. These lines stem from the decay into 7/1 and 717 and would appear at different 
energies as long as the higgs and the 77 masses are not too degenerate. In case C both 
of the decay channels are in general open as long as sin /3 7^ 0, which is the case for 
benchmark point 4 in Tab. [Tj In Fig. |4] we show the corresponding cosmic-ray fluxes. 
Most interestingly the gamma-ray flux exhibits two strong peaks in this case, located at 
270 and 770 GeV. 
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Figure 1: Predictions for case A, benchmark 1, with tdm = 1-7 x lO^'^s (A = 2.9 x 
lO^^GeV). The upper panels show the positron fraction {left) and the total electron + 
positron flux {right) compared with experimental data. Dashed lines show the adopted 
astrophysical background, solid lines are background + dark matter signal (which overlap 
the background in this plot). The lower left panel shows the gamma-ray signal from dark 
matter decay, whereas the lower right panel shows the p/p-ratio: background {dashed 
line) and overall flux {solid lines, again identical with background). 



Case D. This operator, see Eq. (2.7), is particularly interesting since it induces a kinetic 
mixing between the U{1)y of hypercharge and one of the hidden SU{2) gauge bosons. 
As a result two-body decay modes into lepton and quark pairs are allowed, in contrast 
to the other operators. This leads to interesting implications for the electron/positron 
flux that will be discussed shortly below. 

Here we firstly emphasize that again the operator also predicts two-body decay into 
7/1, which could be observable in different parts of the parameter space. The inverse 
decay rate reads, for <^ Ma- 

XI 2H f A \ VlTeV\ V300GeV\^ , , 

and shows that the line could be observed by Fermi LAT for scales of the custodial 
symmetry breaking close to the Grand Unification scale. For these large lifetimes around 
10^* s contributions to the anti-matter channel would be negligible. However, if the line 
lies above around 300 GeV and out of reach of Fermi LAT, shorter lifetimes cannot be 
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Figure 2: Like Fig. [T| but for case A, benchmark 2, with tdm = 1-1 x 10^^ s (A = 
3.7 X 10^^ GeV). The yellow band shows the uncertainty in the anti-proton flux due to 
the propagation model parameters. 



excluded and the anti-matter fluxes can be sizeable. 



3.3 Positron flux 

Here we will briefly discuss the predictions for the anti-matter fluxes concentrating on 
case D, since this operator features two-body decay into fermions pairs due to effective 
kinetic mixing between hidden sector and the hypercharge f/(l)y. The corresponding 
branching ratios are listed in Tab. |4j In the cases with lower dark matter mass, the 
branching ratio into hard leptons (and in particular electrons) is sizable. Namely, in the 
limit Ma ^ Mz the inverse decay rate into charged lepton pairs is given by 

which can produce a steep rise in the observed positron fraction for values of the scale 
of custodial symmetry breaking of the order of the Grand Uniflcation scale. 

As an example, we show in Fig. |5] the predictions for the cosmic ray fluxes for 
benchmark point 2. For a scale of custodial symmetry breaking A = 7.2 x 10^^ GeV, which 
is close to the Grand Unification scale, the gamma ray spectrum shows a intense gamma- 
ray line at 300 GeV, in agreement with current observations. On the other hand, the 
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Benchmark 


Zr\ 


7// 


Zh 


'-yh 


1 


0.19 


0.81 








2 


0.22 


0.78 








3 


0.23 


0.77 








4 


0.028 


0.79 


0.041 


0.14 



Table 3: Branching Ratios for Case C, including benchmark point 4 which features decay 
channels with h in the final state. 



0.050 




Energy [GeV] liin- Energy [GeV] 

Figure 3: Like Fig. [l| but for case C, benchmark 3, with tdm = 6.0 x 10^^ s (A = 
2.0 X lO^^GeV). 

positron fraction shows a steep rise which could partially, although not totally, contribute 
to the PAMELA positron excess. Moreover, the decay into charged leptons is necessarily 
accompanied by a decay into quarks, which produce a sizable antiproton flux and is in 
some tension with the observations. This is a generic feature of the decay mode and hence 
it is unlikely that it contributes the dominant part to the observed positron excess. 

In more generality we found that the PAMELA and Fermi results can be reproduced 
in principle by the model, but only at the price of producing a too large diffuse 7 signal, 
too many antiprotons (unless the dark matter is very heavy) and sometimes gamma lines 
above the rates allowed by the H.E.S.S. measurements in the multi TeV range. 
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0.050 




Energy [GeV] kin. Energy [GeV] 

Figure 4: Like Fig. [l| but for case C, benchmark 4, with tdm = 1-6 x 10^'' s (A = 
1.2 X 10^6 GeV). 
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w+w- 


vv 




uu 


dd 
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0.01 


0.005 


0.04 


0.02 


0.09 


0.39 


0.29 


0.15 
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0.019 


0.004 


0.036 


0.014 


0.072 


0.35 


0.39 


0.12 


3 


0.22 


0.0002 


0.73 


0.0005 


0.003 


0.016 


0.018 


0.005 



Table 4: Branching Ratios for Case D. 



Discussion. It is intriguing that the production of a 7-ray hne is a generic prediction 
for all possible operators that may mediate the decay of the 5'[/(2)hs dark matter gauge 
bosons. For values of the custodial symmetry breaking scale near to the Grand Unification 
scale, and for dark matter masses around 400 GeV and below, this line could be in reach of 
sensitivity of the Fermi LAT gamma-ray line searches. On the other hand, a production of 
an observable amount of electrons and positrons or anti-protons is very model dependent. 
In most cases electrons and positrons are produced in the fragmentation of scalar or 
vector bosons and lead to a very flat spectrum. An interesting exception occurs for 
the operator case D which features two-body decay modes into lepton pairs. In this 
case the produced positron spectrum can rise more steeply, but, when also taking other 
observations into account, still not enough to explain the PAMELA observations alone. 
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Figure 5: Like Fig. [T| but for case D, benchmark 2, with tdm = 6.7 x 10^^ s (A = 
1.5 X 10^^ GeV). 

4 Effects of the annihilation processes with one dark 
matter particle in the final state 

The model considered above has the interesting and rather pecuhar property that it 
allows annihilation processes with one dark matter particle in the final state, i.e. AiAj — )■ 
AkT] annihilations via an intermediate Ak, Fig. 2 of [1]. In ordinary models based on a Z2 
symmetry such processes are strictly forbidden, they would be equivalent to Z2 breaking 
at the renormalizable level and therefore to fast DM decay. The non-abelian character 
of the custodial symmetry responsible for the stability of the hidden vectors allows these 
processes through the trilinear coupling C 3 —^F^^^F^^ 3 —}^gfj,eijkA^A^f^{d^Aiy — dyAi^). 
As pointed out in Refs. [U |2] these "trilinear" processes do not bring nevertheless any 
new radical change in the freeze-out mechanism. In the Boltzmann equations (where 
n = rii + + ^3 is the density of A states) 



dn 



+ 3Hn 



{(^iiV) I 2 



( 



n 



n 



Eq 



) 



{(JijV) 



n[n 



riEq) 



(4.1) 



these terms, parametrized by cTjj, behave linearly in n — n^;^, whereas the ordinary anni- 
hilations, parametrized by an, behave linearly in — n'^Eq- Since — n'^Eq ^ '^n{n — nEQ) 
near freeze-out, the relic abundance behaves as usual ^dm oc l/Max((Tjj, 2(Tjj). However 
these "trilinear" processes contribute with a rate expected to be similar to the one of 
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Figure 6: Parameter space leading to 0.091 < Qh'^ < 0.129 for the small regime 
(10"'^ < Xm < 10~^). From left to right: vs g^, Mr, vs Ma and Ma vs g^. The 
different curves are for various values of A^: Xcf, = 10~^ (red), A,^ = 10"^ (orange), 
A<^ = 10~^ (green) and A^ = 10~^ (blue). The Higgs mass is fixed at Mh = 120 GeV. 
The dots off the main "sequences" correspond to Higgs or rj resonant annihilations, for 
Ma = g(f>V(ii/2 ~ M/i/2 and M^/2 respectively. 

the processes with no dark matter particle in the final states and consequently should be 
properly taken into account. This is what is done here, for "small" Higgs portal coupling, 
Xm < 10~^, and for "large" Higgs portal coupling. Am > 10~^. 

Considering first the small A^ regime, in Fig [6] are shown the values of the gauge 
coupling g^ vs v^, Ma vs and Ma vs g^, which lead to a relic density within the 
WMAP range 0.091 < fi/i^ < 0.129 at 3a |68], in agreement with the direct detection 
experimental upper bounds from CDMS [69] and XenonlO [70] (see Fig. g below). These 
graphs show corrections of order unity with respect to the corresponding result without 
the "trilinear" processes. Fig. 3 of Ref . [1] . The dominant processes are the annihilations 
AiAi — )■ rjrj and AiAj — )■ rjAk which, unless the A<^ coupling is large, have cross-sections 
proportional to ^f^/M^, leading to a Ma oc g'^ quadratic behavior in Fig. [6j The only 
exception to this behavior is given by the resonant cases, when Ma ~ Mh/2 or Ma ~ 
M^/2. 

The corresponding plots for the large Higgs portal regime are given in Fig [7j In this 
case the deviations due to the new AiAj — )■ A^rj processes are more difficult to single out, 
since more annihilation channels (involving Am) contribute to the relic density. But with 
respect to the case already discussed in Fig. 4 of Ref. [1], one finds points with lighter 

and smaller g^ for a same value of f<^. The plot Ma vs g^j) indicates again that the 
freeze-out has a complicated dependence on the couplings of the model. Some of the dots 
still display the quadratic behavior of Fig. |6| when the dominant annihilation channels 
are A^Ai — )■ rjr] and A^Aj — )■ AkT]. 
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Figure 7: Parameter space leading to 0.091 < Qh'^ < 0.129 for the large regime 
(10~^ < Am < 1), from left to right: vs g^, vs and Ma vs g^f,. is varied 
in the range 10~^ — 1, the Higgs mass between Mh = 114.4 GeV and Mh = 160 GeV. 
Here too one can recognize the resonant cases, from both the Higgs and the t] bosons, 
for Ma = g^v^/2 ~ m/j/2,m^/2 respectively. All dots satisfy the LEP constraints on the 
T parameter and on h ^ ff decay. They also are in agreement with the direct detection 
bounds from CDMS [69] and Xenon [70]. The black dots correspond to instances which 
lead to a spin-independent elastic cross-section at most one order of magnitude below 
these bounds. 



At tree level the elastic scattering of the vector dark matter on a nucleon n is spin 
independent, mediated by an /i or boson, and the full expression for the cross-section 
reads: 



1 ,24^2,. 



o-iAn An) ^ ^J%.n.'20 , (4.2) 

with /i„ = ninMA/ (m„ + M^) the reduced mass and m„ the nucleon mass. The parameter 
/ designs the Higgs nucleon coupling, / = {n\ '^qrnqqq\n) and is taken to the value of 
/ = 0.3. 

Imposing the relic density constraint, the predictions for the direct detection rate 
are given in Fig. |8| together with the upper bounds of CDMS [69], XenonlO [TO] and 
the recent final results from CDMS-II [71]. In the small Higgs portal regime. Am < 10^^, 
even though cross sections are suppressed by 2 powers of A^, large direct detection rates 
can be obtained for small A,^ couplings because in this case < m/j and the An — )■ An 
cross section scales as A,^/A^, Eq. (4.2). For a dark matter mass from few GeV all the 
way up to the multi TeV range, a spin independent cross section of the order of the 
current experimental sensitivities can be obtained for values of A<^ of order or below few 
10"'^. The values of the various parameters required in this case can be read off in Fig. [6] 
The value of the cross-section does not depend on the dark matter mass if this mass is 
large, as indicated by the plateaux for different values of A^. For the large Higgs portal 
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Figure 8: Obtained spin independent cross-section on nucleon a^^{An An) versus M4, 
in agreement with the constraint 0.091 < Qh"^ < 0.129. Small regime (10^^ < < 
10^'^) on the left and large Higgs coupling portal, A^ > 10~^ on the right. The color 
caption is as in Figs. [6] and [7j The thick (dashed) black curve is the CDMS (XenonlO) 
upper bounds at 90% C.L.. The dotted-dashed curve is the recent published CDMS-II 
upper bound, at 90% C.L. 

regime. Am ^ 10^^, direct detection rates of order the present experimental sensitivity 
or exceeding it are easily produced. For illustration, among the sets of parameters that 
lead to the right relic density in Fig. [7| we have denoted by black dots the ones which 
lead to an elastic cross-section on nucleon at most one order of magnitude below the 
CDMS [69] and Xenon [70] limits. Here too a dark matter mass in the whole range from 
1 GeV up to few TeV can be accommodated. Even though other values are possible, the 
rj mass tends to be either small, below 100 GeV, or slightly larger than the dark matter 
mass. For larger than ~ 700 GeV one recovers the linear relation between and 
and the corresponding quadratic behavior of Ma in 5^^, indicating that the pure hidden 
sector annihilations driven by the coupling are dominant, as in the small Higgs portal 
regime. 

Note that all the dots shown in the figures above satisfy the LEP constraints. The 
mixing of the rj boson with the standard model Higgs affects the electroweak precision 
observables. The main constraint on the model parameters comes from the T parameter, 
since the 77 is a neutral scalar which mixes with the Higgs boson. We use the same cuts 
as in Ref. [Ij, that is to require that T — Tsm is in the conservative range —0.27 =F 0.05 
from [72]. For < 114.4 GeV the branching ratio rj ^ ff should not exceed the LEP 
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direct search bounds, leading to an upper value on the mixing angle sin^ /5, see Fig. 10 
of Ref [73]. This bound constraints more the large Higgs portal coupling regime which 
involve large mixing angles. In addition, we take into account the combined analysis of 
the CDF and DO collaborations, that excludes a Higgs in the mass range 160 — 170 GeV 
at 95% C.L. [71]. 

5 Conclusions 

We have shown in this paper that the vectors of a hidden, spontaneously broken, non- 
abelian gauge group constitute a viable dark matter particle which decays at cosmological 
times. Their longevity is due to an accidental custodial symmetry in the renormalizable 
Lagrangian. However, similarly to the proton, they are not expected to be absolutely 
stable due to the existence of non-renormalizable dimension six operators which induce 
the decay of the dark matter particle. We have identified the four dimension six operators 
which break the custodial symmetry and calculated the dominant decay modes. Taking 
advantage of the fact that the dark matter has spin-1, the gamma lines are produced at 
tree level from DM two-body decay which beside the 7 involve the standard model h or 
hidden sector r] Higgs bosons. We have found that in all the cases an intense gamma-ray 
line is expected, which could be observed by the Fermi-LAT if the scale of custodial 
symmetry breaking is close to the Grand Unification scale. We have also calculated the 
positron fraction, total electron plus positron flux and the antiproton-to-proton fraction 
for these channels. Even though in these scenarios there is a sizable branching ratio into 
hadrons, the total antiproton-to-proton fraction is consistent with the measurements, 
while still producing an observable and possibly huge gamma-ray line. Finally we also 
improved the calculation of the relic density of hidden vectors including the dark matter 
annihilation processes with one dark matter particle in the final state. Direct detection 
rates close to the present limits can easily be produced for any dark matter mass within 
the GeV-multi TeV range. 
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